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The grade of service of a -probability-engineered trunk group is defined 
to be the 20-day average blocking during the busy hour of the busy season. 
In this paper, an improved model for calculating grade of service is de- 
veloped and used to increase the accuracy of the existing trunk-engineering 
procedures. Using the new model, new traffic-capacity tables and trunk- 
estimation algorithms have been designed for use in the Bell System. 

I. INTRODUCTION 

The grade of service for a probability-engineered trunk group is 
defined to be the average blocking observed in the time-consistent 
busy hour of the busy season.* The existing methods for predicting 
grade of service do not account for the effects of the finite length of the 
individual one-hour measurement intervals and thereby tend to under- 
estimate trunk-group capacity. 

In this paper, we develop an improved model for calculating average 
blocking that includes the two essential effects of the finite measure- 
ment interval. First, the current method for estimating the mean 
blocking for a single hour must be revised to remove an implicit 
assumption that the measurement interval is infinite. Second, the 
existing mathematical model for day-to-day variation of trunk-group 

'Both the CCITT (International Telegraph and Telephone Consultative Com- 
mittee) and the Bell System define grade of service as an unweighted average of 
busy-hour busy-season blocking values. 
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offered loads must be modified to account for statistical measurement 
error, which is also introduced by the finite measurement interval. 

We develop a new model of day-to-day load variation in Ap- 
pendix A and in Section II combine it with a new estimate of mean 
single-hour blocking to obtain our approximation for the average 
blocking. This new approximation is then compared with the existing 
approximation analytically and numerically; the accuracy of the new 
approach is established in the third section using data from a computer 
simulation. A summary is given in the last section. 

II. AVERAGE BLOCKING 

In this section, we develop an approximation for the average busy- 
season busy-hour blocking. Since the measured grade of service is de- 
termined from blocking measurements made over several hours, and 
since the busy-hour source loads vary from day to day, our analysis 
must account for the effects of such load variation. Accordingly, we 
first discuss day-to-day load variation and then describe our approxi- 
mation for average blocking. 

2.1 Day-to-day load variation 

R. I. Wilkinson was the first author to study the impact of day-to- 
day variation in offered loads on trunk-engineering procedures. 1,2 He 
collected data from a number of trunk groups that indicated that the 
distribution of the observed loads could be approximated by a gamma 
distribution. The data also indicated that the variance, Var (&), of 
the observed load, &, was related to the mean a by 

Var (&) = 0.13a*, (1) 

where is a parameter whose value depends on local conditions. 

Wilkinson's studies showed that Var (d) tends to be relatively larger 
for overflow traffic than for first-routed (Poisson) traffic. Those results 
led to the specification of three values of <f> (1.5, 1.7, and 1.84) to cover 
a reasonable range of engineering applications. The level of day-to-day 
variation is called low when <j> = 1.5 is appropriate, medium for 
<f> = 1.7, and high when <f> = 1.84. 

In Appendix A, we show that the variance of the observed "single- 
hour" offered loads is a sum of two components : 

Var (&) = -^ + Var (a), (2) 

where a is the daily source load (a random variable), & is the observed 
load, 5 = E{a) is the average busy-season load, z is the traffic peaked- 
ness, t is the observation interval (usually one hour), and h is the 
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mean holding time. The first term is an approximation for the variance 
arising from the finite measurement interval* and Var (a) is the vari- 
ance of the daily source load. 

Combining eqs. (1) and (2), we obtain a model relating the source- 
load variation to the observed load variation; i.e., 

Var (a) = 0.13a* - ^EL- (3) 

For certain combinations of the various parameters, 2az/ (t/h) can ex- 
ceed 0.13a*, indicating that the observed load variation is entirely 
due to random measurement error resulting from the finite measure- 
ment interval. For our application, we will assume that 



Var (a) = max 



°.°- 13a *-<i))- < 4 > 



2.2 The approximation 

Consider a service system with c servers having exponentially dis- 
tributed service times with mean h and serving traffic under a blocked- 
calls-cleared service discipline. The system is observed during n dis- 
joint measurement intervals 1%, ■ ■ -, I n , each of length t. During /*, 
the interarrival times are independent and identically distributed (iid) 
with mean 1/A*. The peakedness, z, of the traffic is assumed to be the 
same during all the intervals. + The system is in statistical equilibrium 
during each interval and the initial point of each interval is a stationary 
(random) point for the arrival process. The loads a, = X,7t, i = 1, • • • , n, 
are independent and identically distributed according to the distribu- 
tion function r(a|a, Vd) with mean a and variance Vd = Var (a) (the 
day-to-day source-load variance). [We assume that F(a\a,Vd) is a 
gamma distribution. Justification for the assumption is given below.] 

We use Aj(t) and Oj(t) to denote, respectively, the number of arrivals 
(call attempts) and the number of blocked attempts (the overflows) 
during /,■; the (measured) observed blocking is B, = Oj(t)/Aj(t) pro- 
vided Aj(t) 9± 0. If Aj(t) = 0, then Oj{t) = and the ratio is not 
defined. However, if no arrivals occur during /,-, it seems appropriate 
to say that no blocking occurred; i.e., we define Bj = whenever 
Aj(t) = 0. The sample average of the observed blocking (the observed 
grade of service) is 

1 n 
B n = - Z Bf. 

7lj = i 



* This is an extension of the formula given by Riordan in Ref. 3. 

* Traffic studies have shown that the peakedness of the traffic offered to a final 
trunk group does not change significantly from day to day during the busy hour of 
the busy season. 
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Since a lf ••-,«„ are iid, B lt ■ • •, B n are also iid. Consequently, the 
average blocking B is given by 



B = E(B n ) = E{B,). 



(5) 



The hourly measurement Bi represents a random sample of observed 
blocking corresponding to a load population with mean ct\. The 
parameter a x is a random variable with mean a and distributed accord- 
ing to T(ai\a, v d ). Thus, 



/} 



"/. 



Oi(t) 



9 E 1 Al(0 



0L\ 



dT{ai\a, v d ). 



(6) 



Dropping the subscripts, the conditional mean in (6) is given by 



E 



0(0 



AM 



= Pr {A(t) > 0\a\E 



0(0 



A(t) 



a, A(t) > 



where Pr (X) denotes the probability of the event X. Let 
A = E{A(t)\a\ and = E{0{t)\a). Also, for nonnegative m and n, 



^{[O(0]-[ii(0]-|«,A(0>0} = 



E{lOm m lA{t)y\a] 
Pr [Ait) > 0|a} 



Using these relations, an approximation for the conditional expecta- 
tion is obtained by expanding the function x/y in a two-dimensional 
Taylor Series about the point (x , yo) = (E{0(t)\a, A{t) > 0}, 
E{A(t) \a, A(t) > 0}). Taking the appropriate conditional expectation 
and retaining only the terms up through second order, we have the 
approximation 



E 



0(0 



A{t) 



OL, A(t) > 



. O/Pr {A(0 > 0|a} V 



5 + 2( 

AAV 






r E[\jA{t)j\ a ) ( 

|_Pr \A{t) > 0\ a ] \Pt {A(t) > 
_ / Pr jA(t) > 0\a] \ 2 

/ E{A(t)Q{t)\a) AO \ 

' \Pr {A(t) > 0\a\ (Pr {A(t) > 0|a}) 2 / 



(7) 



Denoting the call congestion O/A by B(c, a, z) and noting that 
E{A(t)\a) = at/h, eq. (7) reduces to 



E 



O(0 



A{t) 



a,A(t) > ttB(c, a, z) + 



Pr \A{t) > Ola} 
(at/h) 2 

lB{c, a,z) Var {A(0|a} - Cov {A(t),0(t)\a}2, 
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where all of the statistical moments are functions of a, z, and c. Thus, 
Oil) 



E 



A(t) 



a\ ^Pr \A(t) > 0\a}B(c,a,z) + ( Pf &&> 9M )' 
•ZB(c,a,z)Var{A(t)\a} - Cov {A (*), (<)!«}]. (8) 



Formulas for computing the moments are given in Ref. 4. An expres- 
sion for Pr {A(t) > 0\a] is given in Appendix B. Numerical experi- 
mentation has shown that Pr [A(t) > 0\a\ ~ 1 for at/h > 10; i.e., 
that term can be ignored except for very small loads or short measure- 
ment intervals. Combining eqs. (6) and (8), we have 

Btt ["PrUW > 0\a\B(c, a, z)dT(a\a, v d ) 
Jo 

+ ^ R(c,a,z)dT(a\a,v d ), (9) 
where 

• [£(c, a, z) Var \A{t)\a] - Cov {A(t), 0(t)\a}J 

The approximation is complete when V is specified. 

Numerical experimentation has shown that B is not very sensitive 
to the shape of T(a\a, v d ) for fixed values of 5 and v*. Accordingly, 
following Refs. 1 and 2, we have assumed that T(a\a, Vd) is a gamma 
distribution with mean a and variance v d , where v* = Var (a) is given 
by eq. (4). With these assumptions, the integrals in (9) can be com- 
puted numerically using a 51-point compound Simpson's rule. The 
accuracy of the approximation is discussed in Section III. In Section 
2.3, we compare this approximation with the existing procedure. 

2.3 The existing procedure 

The existing approximation for average blocking is 1,2 

B e (c, 5, z) ^ f° B(c, a, z)dT(a\a, v), (10) 

Jo 

where 

v = Var {&) = 0.135*. 

Consequently, there are two essential differences between the present 
method (10) and our approximation (9). 

First, the existing method neglects the random component of load 
variation and assumes that all of the variation in the observed loads 
is due to day-to-day changes in the source load; i.e., the assumed load 
variation is too large by an amount 2az/(t/h). The integral in (10) is 
an increasing function of v in the range of engineering interest (less 
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than 5 percent average blocking) and, hence, B e is larger than it would 
be if va were used in place of v. 

Comparing (9) and (10), we see the second difference. The existing 
procedure implicitly uses E[[0(t)/A(t)2\a} Prf B(c, a, z) and neglects 
the term R(c, a, z). Numerical experimentation has shown that 
R (c, a, z) is negative in the range of engineering interest and that 
\R(c,a,z)\ is an increasing function of z. For z near 1, R(c, a, z) 
is negligible. However, for z ^ 2 R(c, a, z) becomes significant. Con- 
sequently, B(c, a, z) is larger than E{[0(t)/A(t)2\a\, and the differ- 
ence increases with z. 

The two differences combine to cause B < B e in all regions of 
engineering interest; i.e., the average of the observed single-hour 
blocking probabilities is less than that predicted by the existing method 
(10). Quantitative comparisons are made in Section III. 

III. NUMERICAL RESULTS 

To determine the accuracy of our approximation (6) relative to the 
existing procedure (10), a computer simulation was constructed for a 
full-access trunk group satisfying the assumptions specified in Section 
2.2. The simulation was run for a large range of the system parameters 
c, a, z, and <f> covering the regions of engineering interest (B tt 0.01, 
c ^ 2, 1 ^ z ^ 7). Generally, a mean holding time of 180 seconds 
was used (although both smaller and larger values were used for sen- 
sitivity tests)* and 20-day averages were generated. All statistics are 
based on a sample size of 50; i.e., 1000 simulated hours. 

3.1 Simulation output 

Typical results from the simulation are summarized in Tables I and 
II. First consider Table I. The first four columns are the input parame- 
ters for the simulation. In order, they are the peakedness z, the offered 
load a, the trunk-group size c, and the conventional exponent <j> defining 
the level of day-to-day variation. QThat is, Var (a) was adjusted so 
as to produce the desired <f>, as discussed in Section 2.2.] The next 
three columns give the simulated 20-day average blocking B , the 
variance v a of the daily offered source loads, and the (total) variance 
v a of the observed loads. The measurement variance vs is then com- 
puted as Va = v- a — v a , and is compared with 2az/(t/h) in the last two 
columns. In the cases shown, 2az/(t/h) is an adequate approximation 
of the variance introduced by a finite measurement interval. The other 



* The holding time was varied from 100 to 360 seconds and a small effect was 
observed. The effect is negligible for most engineering applications. Based on a survey 
of observed holding times, AT&T has requested that h = 225 seconds be used for 
the new traffic tables. 
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Table I — Simulation results 



Simulation Input 


Observed Data 


Peaked- 


Mean 


Trunks 
(c) 


Variation 


Average 


Input 


Observed 


Residual 


Theoretical 


ness 


Load 
(5) 


Parameter 
(*) 


Blocking 
(B ) 


Variance 
(v a ) 


Variance 

(va) 


Variance 
(i>5) 


m 


1.0 


4.01 


10 


1.5 


0.0083 


0.72 


1.03 


0.31 


0.40 


4.0 


17.80 


40 


1.5 


0.0084 


4.81 


12.54 


7.73 


7.12 


4.0 


9.80 


30 


1.7 


0.0049 


3.60 


6.61 


3.01 


3.92 


7.0 


9.75 


40 


1.84 


0.0046 


1.76 


8.96 


7.10 


6.82 



test cases indicated that the accuracy of the approximation generally 
increases as a increases (as one would expect from the asymptotic 
nature of the approximation). 

3.2 Existing approximation 

The data in the first seven columns of Table II illustrate the size 
of the bias in the existing method. The first five columns of the table 
are the same as those in Table I. The next column B e represents the 
computed average blocking corresponding to the source load, peaked- 
ness, trunk-group size, and respective day-to-day variance used in the 
simulation. The next column c e gives the corresponding estimate (using 
the existing method) of the trunk-group size necessary to achieve the 
simulated blocking for the given input load, peakedness, and day-to- 
day variation parameter <f>. 

First, note the difference between B and B e . The existing approxi- 
mation B e is always larger than the actual average blocking B , and 
the relative difference increases as z increases. The bias in B e will 
cause the engineering estimates of trunks required (to meet objective 

Table II — Engineering methods 



Simulation 


Engineering Methods 


Input 


Output 


Existing Method 


New Method 


Peaked- 
ness 

w 


Mean 

Load 

(«) 


Trunks 
(c) 


Variation 
Parameter 

(*) 


Simulated 
Average 
Blocking 

(So) 


Blocking 

Estimate 

(B e ) 


Trunk 

Estimate 

(c.) 


Blocking 

Estimate 
(B) 


Trunk 

Estimate 
(c) 


1.0 
4.0 
4.0 
7.0 


4.01 

17.80 
9.80 
9.75 


10 
40 
30 
40 


1.5 
1.5 
1.7 
1.84 


0.0083 
0.0084 
0.0049 
0.0046 


0.0100 
0.0145 
0.0103 
0.0106 


10.24 
42.83 
32.96 
44.71 


0.0078 
0.0083 
0.0054 
0.0045 


9.91 
39.97 
30.31 
39.93 
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service) to be too large. The bias in the trunk estimates is illustrated 
by comparing c with c e .* 

In all the cases we considered, the bias was primarily a function of 
z. For 2=1, the bias was generally less than one trunk. For 2 = 2, 
c e — c was usually between one and two trunks, and, for 2 = 4, the 
bias ranged between throe and five trunks. The corresponding relative 
errors (c — c e )/c were largest at the smaller values of c. 

3.3 New approximation 

The relative accuracy of the new approximation is illustrated in the 
last two columns of Table II. In these cases, the new approximation 
for B is much closer to the simulated blocking B . The relative differ- 
ences are quite small, especially when compared with the correspond- 
ing errors in the existing approximation. Since the estimate of B is 
good, the corresponding estimate c of trunks required to achieve B 
is quite close to c, the number actually required. 

Similar results were obtained in all the other test cases. Accordingly, 
we conclude that the existing approximation (10) is biased, but the 
bias is essentially removed by using the new approximation (9). 

Based on these results, the new approximation has been used to 
generate new trunk-engineering tables and algorithms for use by the 
Bell System. 

IV. SUMMARY AND CONCLUSIONS 

By combining a new model for day-to-day load variation with a 
new estimate of mean single-hour blocking, a new approximation was 
obtained for estimating the grade of service for probability-engineered 
trunk groups. Using this result, it has been possible to improve the 
accuracy of the presently recommended trunk-engineering procedures 
and thereby realize an increase in predicted trunk-group capacities. 
The increases are smaller for trunk groups serving Poisson traffic, 
but become substantial as the peakedness and level of day-to-day 
variation increase. The new approximation has been used to develop 
new trunk-engineering tables and algorithms that will soon be intro- 
duced into the Bell System. 

APPENDIX A 
Day-to-Day Load Variation 

Let S denote a GI/M/<x> service system which is observed over k 
disjoint intervals of time \I\ } • • •, I k ), each of length t.* The inter- 



* Refer to Section 2.3 for a discussion of the individual components of the difference. 

' The infinite server system is a standard model for traffic engineering. It is used to 

characterize a traffic process in terms of only interarrival times and service times. 
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arrival times during /,■ are independent and identically distributed 
(iid) according to the probability distribution function F, with mean 



h-l" urM - 



The arrival rates {Xi, • • •, X*} are iid with mean E{\) and variance 
Var (X). The peakedness (variance-to-mean ratio) of the arrival process 
is constant for all intervals. The service times are iid according to a 
negative-exponential distribution with mean h. The system is assumed 
to be operating in statistical equilibrium at the beginning of each 
interval; i.e., the initial point of each interval is a stationary point for 
the arrival process (see Ref. 3). 

Let Nj be the number of arrivals during /,■ and let £,•> denote the 
service time of the tth arrival in /,-. The average usage during 
/ = 5Z*«i /,• is estimated by* 

1 k Ni 

d* = £ E E fry. (11) 

K y=i ,=i 

(We assume the edge effects are negligible since S is in equilibrium 
during each interval //.) The corresponding estimate of average offered 
load is 

dk-jth. (12) 

In this section, we obtain the mean and variance of &k. 
From Ref. 5 and eqs. (11) and (12), we have 



*™-h&%[ E %M N '}\ 



= - L E [hNi\ 



"B £*«*'»■ (13) 

Since the beginning of each interval is a stationary point for the arrival 
process, it follows that E{Nj\\ } } = \jt (see Ref. 6). Thus, for the 
arrival process, it follows that 

E{N,} = E{E{N;\\j}} 

Xy 

= E [\jt] 

x, 

= E{\}t 



* Our studies have shown that the additional measurement variance caused by 
discretely sampling the usage with a 100-second-scan Traffic Usage Recorder was 
negligible when compared with the variance caused by the load variation. 
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and so from (13) we have 

E{& k ) = E{\}h. 

Since the fc> are independent, eqs. (11) and (12) yield 



Var(<3*) = ^EVar 



E ta\ , 

t = l 



(14) 



and from Ref. 5, 

v -|£M-?te'U b ' i m + i"feUH)- (l5) 



We first expand 



Var< 

Ni 



Ni 

E fc/IX* 

t=i 



= VarM j? foli^X,}} 

+ |{var{ , E fc,|iVy,X,}} 
= Var {iV^IX,} + JS {tfjA*|Xy} 

AT; AT/ 

ui ,/i , VarJiVylX^X 



(16) 



Let 2 denote the traffic peakedness (which is constant over I). Then 
from Ref. 7, we have the approximation 

Var {JVy|Xy} 



22 - 1 ; 



(17) 



B[N,\\,) 

the approximation has been found to be quite good for a > (2 — 1) 
and t ^ lOh (and probably acceptable for engineering purposes for 
t ^ bh). Substituting (17) into (16), we obtain 



Var 



Ni 



E b#|X,[ = 2ehH\i. 



(18) 



Expanding the second term in (15) in the same manner as in eq. 
(13) yields 

E €«|X/} = ht\i. (19) 



E < 



Substituting (18) and (19) in (15) provides 

= 2zhHE{\\ + (fa) 8 Var {X}, 
which, with eq. (14), implies that 



Var j E in 

i = l 



(20) 



Var(^)=^(^^ + / i 2 Var{X})- (21) 
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The offered load during Ij is ay = Xjh. Hence, E{aj\ = hE{\] and 
Var {a,} = h 2 Var {X}. Dropping the subscript./, we have 

v «f*>=K ? f!r L+Var,a1 ) - (22) 

Equation (22) has a simple interpretation. The first term (2zE[a\)/ 
(t/h) represents the component (of the observed variation) which is 
due to a finite measurement period. The second component Var {a} 
is the "true" (day-to-day) variation of the offered source load. The 
factor k results from averaging k observations. Setting k = 1 and 
a. = E{a}, we have the variance of the single-hour load estimate 

Var {a} = || + Var {a}, (23) 



,vhere the individual terms have the interpretations noted above. 



APPENDIX B 

Probability of an Arrival 

In this appendix, we derive an expression for Pr [A(t) > 0|a}. We 
assume that the interarrival times are independent and identically 
distributed according to the distribution function F. We further assume 
that F is a mixture of exponentials, with the parameters chosen so 
that F approximates the interarrival distribution of an overflow pro- 
cess with load a and peakedness z (see Ref. 8) ; that is, 

F(t) = 1 - kxe-^ 1 - k 2 e- rit . 

Now, Pr [A(t) > 0\a\ is just the probability that the first arrival 
after a random entry point (i.e., a stationary point for the arrival 
process) will occur before t units of time have elapsed. Thus, from Ref. 
9, we have 

Pr \A{t) > 0\a\ = a /"' [1 - F(x)~\dx 

Jo 

\ ri r 2 ri r 2 / 

Since f x [1 - F(x)]dx = l/a, it follows that 

Pr [A(t) > 0|a} = 1 - « (y e-"' + ^e-'A , 

where the parameters ki and r» are functions of a and z as specified in 
Ref. 8. 
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